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ABSTRACT 


A method  for  determining  an  optimum  safe-life  is  presented. 

This  life  is  optimal  in  the  tense  that  the  combined  expected  costs 
associated  with  catastrophic  failure  and  total  life  inventory  are 
a minimum. 

An  associated  question  that  is  addressed  is  that  of  determining 
the  optimum  number  of  tubes  to  be  life-tested.  This  decision  is 
resolved  by  minimizing  the  combined  cost  of  testing  and  the  monetary 
penalty  of  imperfect  information. 
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INTRODUCTION 


The  question  of  how  long  to  retain  gun  tubes  in  service  (in 
terms  of  cumulative  usage)  has  historically  been  settled  by  deter- 
mining the  number  of  rounds  at  a givfen  zone  that  a gun  would  fire 
with  acceptable  accuracy.  Due  to  several  factors , nearly  all  the 
cannon  tubes  in  our  inventory  now  become  dangerously  liable  to  a 
fatigue  failure  before  wear  causes  an  unacceptable  loss  of  accuracy. 

Thus,  the  question  of  current  concern  is  how  many  rounds  may 
be  fired  before  the  risk  of  a fatigue  failure  becomes  excessive. 
Unfortunately  the  variability  associated  with  failure  by  fatigue  is 
quite  large.  With  no  precise,  non-destructive  way  of  ascertaining 
when  e tube  will  fell,  one  only  has  recourse  to  determining  tube 
life  an  the  basis  of  failure  date  gethered  on  the  tube.  Further  the 
large  cost  of  tubes,  used  expeodebly  in  testing,  precludes  large 
sample  tlzea.  Thus  the  risk  that  one  is  running  must  Itself  be 
obtained  by  statistical  inference. 
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Stated  more  precisely,  from  a small  sample  of  failure  data,  one 
must  decide  what  the  nature  of  the  failure  distribution  is  and  what 
the  parameters  of  this  distribution*  are.  Only  when  the  latter  are 
obtained,  can  one  predict  the  chances  of  failure  when  operating  a tube 
to  a certain  life. 

Presently,  statistical  data  on  thick-walled  cannon  tubes  indicates 
that  the  distribution  of  cycles  to  failure  is  log-normal.  Altho  some 
other  distributions  cannot  be  definitively  rejected  by  chi-squared  and 
Kolaogorov-Smirnov  tests,  graphical  analysis  using  maximum  likelihood 
estimates  offers  persuasive  evidence  that  the  distribution  is  log-normal 
rather  than,  say,  Weibull.  Host  of  the  failure  data  from  which  this  in- 
ference is  drawn  cocoa  from  one  tube— the  M113  cannon  used  an  the  M107, 
175mm  gun.  However,  because  of  the  similarity  In  geometry  of  tubes  with 
changes  In  caliber,  It  is  presently  felt  that  the  distribution  of  tube 


* Even  to  talk  about  a failure  distribution  implies  immutability  of 
expected  values  over  time.  This  in  turn  requires  that  the  process 
generating  the  tubes  does  not  change  over  time.  For  example  the 
ssste  materials  must  be  used  end  the  same  requirements  for  strength, 
hardness,  toughness,  etc.  must  persist.  The  same  fraction  of  tubes 
must  be  procured  from  each  supplic..  ,'hi  seme  inspection  standards 
must  apply,  etc.  Doubtless  pver  the  life  of  the  systen  aooa  charac- 
teristics of  the  process  will  change  with  possible  change  in  the 
nature  of  the  product.  This  prospect  adds  another  element  of  uncer- 
tainty to  an  already  incertain  situation.  However,  we  shall  assume 
statistical  statlonarlty,  mainly  because  the  problem  would  become 
completely  intractable  without  this  assumption,  but  also  because  we 
feel  that  changes  deliberately  introduced  in  the  process  would  only 
tend  to  Improve  the  quality  of  the  product  ttd  therefore  that  the 
fictitious  failure  distribution  obtained  from  early  sampling  would 
be  a conservative  estimate  of  the  present  quality  of  the  product. 
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life  would  be  log-normal  for  all  cannons.  In  this  report  we  assume 
that  tube  life 'is  log-normally  distributed.  The  failure  data  for  the 
M113  cannon  and  graphical  evidence  for  the  log-normality  of  the  under- 
lying distribution  are  presented  in* the  appendix. 

Therefore,  only  the  parameters  of  this  distribution  need  be 
estimated  from  early  failure  data.  The  primary  purpose  of  these  tests 
is  viewed  as  particularizing  the  distribution  of  tube  life.  This 
distribution  will  be  celled  F(t)  and  its  estimate  will  be  F(t)  . 
After  this  has  been  done  e decision  must  be  made,  in  the  light  of 
F(t)  , es  to  what  is  a reasonable  safe-life.  One  way  of  doing  this 
is  to  decide  what  a maximum  acceptable  risk  of  failure  is,  say,  a j 
and  find  the  value  of  t such  that  F(t)  * o . 

Altho  this  suggestion  bee  the  merit  of  conceptual  simplicity,  it 
fails  to  address  how  one  decides  upon  e max  acceptable  risk.  Since 
determination  of  the  letter  is  the  heart  of  the  problem,  it  is  not 
helpful  to  answer  the  question  of  tube  safe-life  by  deferring  the 
queation  of  risk  to  executive  fiat. 
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DETERMINING  OPTIMUM  SAFE-LIFE 


The  author*  believe  that  a clearer  picture  of  the  consequences 
of  selecting  a certain  life  as  "safe"  is  obtained  by  examining  the 
future  costs  associated  with  a selection— the  only  costs  germane  to 
this  decision*  These  costs  are  (1)  inventory  costs  and  (2)  expected 
costs  imputed  to  catastrophic  failures.  The  first  of  these  decreases 
as  tube  safe-life  is  increased  whereas  the  second  Increases  with  tube 
life.  This  situation  suggests  that  the  sum  of  these  costs  may  have  a 
minimum  for  a particular  life.  This  life  is,  then,  regarded  as  optlmel. 
To  find  it,  we  proceed  as  follows. 

Let  c • the  average  coat  of  acquiring  and  maintaining  a tube  in 
inventory 

(A) (c)“  the  average  cost  of  an  operational  failure  of  a tube 

including  all  costs  imputed  thereto  such  as  liabilities, 
cost  of  investigation,  value  loss  of  operational  oppor- 
tunity,  value  loss  of  troop  confidence,  value  loss  of 
military  prestige,  etc. 

NOTE:  A is  the  ratio  of  cost  of  failure  to  tube  cost. 

u * the  expected  total  utilisation  of  the  system  over  its 
operational  Ufa  in  number  of  rounds 
t^  - condemnation  limit  (safe-life)  for  a tube 
F(t)  - failure  distribution  for  tube  life 
Then  the  total,  undiscounted,  expected  future  cost  affected  by  a decision 
on  tL  is  given  by 

total  variable  cost  • (A  c u P(t^)  + c u]/t^  (1) 
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We  call 


S 


0(tL)  ■ variable  cost  / (c  u)  , 

(2) 

hence 

U(tL)  - [A  F(tL)  + l}/tL 

(3) 

It  is  clear  that  minimization  of  0(t  ) with  respect  to  t is 

L L 

equivalent  to  minimization  of  the  total  variable  cost.  In  graphs  1 

and  2 we  have  plotted  U(t.)  vs  t.  for  two  values  of  the  parameter 

A . For  this  purpose  we  chose  F(t  ) as  log-normal  with  parameters 

L 

p * 8.3065  and  a - 0.4425,  this  being  representative  of  the  M113  gun 
tube. 

f*  -tJ^2 

F(t^)  • 1 I e dt  , where 
In  ti  - y 

1 " a (4) 

In  graph  3 is  shown  the  locus  ox  optimum  t^  . In  graph  4 we 

have  cross-plotted  the  optimum  value  of  t^  , t*  , vs  the  parameter 

* 

A . Tha  unreliability,  F(t^)  * i*  *l*o  shown.  It  is  clear  from 

this  graph  that  the  optimal  value  of  tuba  safe-life  becomes  progres- 
sively less  sensitive  to  increments  in  A as  A increases. 

For  the  case  shown,  a tube  cost  of  $25,000  has  been  given  while 
tha  total  imputed  cost  of  a catastrophic  failure  of  $2.5  million  is 
not  unreasonable.  Thus  a reasonable  value  of  A is  100.  At  this 
value  the  tube  safe-life  is  1100  cycles  and  is  quite  insensitive  to 
changes  in  A . It  should  be  noted  that  this  optimal  safe  life  is 


considerably  above  the  400  cycles  currently  in  use  for  the  K113  cannon 
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PROBLEM  OF  TEST  SAMPLE  SIZE 


S'  - 


The  executive  will  choose  a tube  safe-life  o»  the  basis  of  the 
best  indorsation  he  has  available  after  testing.  To  use  the  method 
described  above,  ha  must  set  as  if  his  estimate  of  the  failure  dls- 

A 

tributicn  F(tjP  is  the  one  which  actually  obtains.  He  would  then 
proceed  to  Obtain  bis  estimated  optimal  safe  life  tL  * 

However,  because  of  inevitable  sampling  error,  the  value  of  his 
"optimum"  will  not  in  general  be  the  correct  one.  That  Is,  in  general 
tj  / . Because  of  this  the  total  variable  cost  of  subsequent 

operation  will  be  higher  than  it  might  have  been  had  he  chosen  t*  . 

Letting  U*  * U(t£)  , the  monetary  penalty  associated  with  choosing 

A# 

t^  is  given  by 

penalty  • c u (U(t£)  - 0*)  (5) 

This  is  the  loss  Incurred  because  of  Imperfect  information. 

How  one  esn  reduce  the  expected  difference  in  t£  end  t£  by  test- 
ing more  units.  The  question  as  to  how  many  units  one  should  test  can 
then  be  posed,  for  it  it  clesr  that  one  cannot  go  on  testing  indefinitely 
in  in  attest  to  reduce  hie  uncertainty  concerning  F(t^)  . 

We  answer  this  question  by  insisting  that  the  rational  decision- 
maker will  want  to  minimize  the  expected  aggregate  cost  of  testing  and 
penalty.  This  total  cost  for  n units  in  a total  life  test  is 

total  cost  • penalty  + B c n , (6) 


where  B is  the  ratio  of  cost  of  testing  one  tube  to  the  cost  of  a tube 


la  Inventory,  c . This  value  of  B Is  about  12  for  existing  systems 
and  methods. 

He  define 

V(n)  ■ total  cost  / c . (7) 

Thus 


V(n)  « u G(n,t)  + B n, 

where 


G(n,t)  - U(t)  - U*  . (3) 

It  is  clear  that  the  expected  value  of  7(n)  , which  we  shall 

minimise,  is  given  by 

V(n)  ■ u G(n)  + ! a (9) 

To  find  the  value  of  n for  which  7(a)  la  a minimum,  we  proceed 
as  follows. 

Let  n*  be  such  that 
V(n*)  < V(n*  ♦ 1) 

< V(a*  - 1)  (10) 

®ut»  ■*"  » «•  * e 

V(n  + l)=su  G(n*  + 1)  + B (n  +1) 

or  V(n*  + l)  j u[  G(a*  + 1)  - C(n#)  ) + B + V(n*)  . (11) 

Similarly, 


7(a*  - 1)  • u [ ^(n*  - 1)  - G(a*)  J - B + v(n*) 
From  (10),  (11),  end  (12),  we  can  write 


G(n*  - l)  - cf(n*)  » — » C(n*)  * G(n*  + 1) 

u 


(12) 


(13) 


11 


The  last  expression  is  interesting  in  that  it  shows  that  the 
optimum  sample  size  depends  only  on  the  ratio  of  B to  u , i.e. , 
only  on  the  ratio  of  non-dimensional  test  cost  to  utilization.  Further 
at  the  optimum  n , the  forward  difference  in  G(n)  straddles  B/u  . 

To  determine  G(n)  , the  authors  found  it  necessary  to  resort  to 
MONTE  CARLO  techniques.  In  doing  this,  we  attempted  to  duplicate  as 
closely  as  possible  the  conditions  facing  the  tester. 

£-mpla  of  n failure  times  are  drawn  from  a log-normal  distri- 
bution whose  parameters  are  unknown  to  the  tester.  From  these 
tx  (1  • 1 . ..n),  he  uses  maximum  likelihood  to  obtain  parameter 

A A 

estiuatec  y and  o 

y • En  In  \t.)  / n 
i-l  1 

o - sqrt  [ (I*  In2  (tt)  - nu2)/(n-l)J  (14) 

i-l 

A 

The  estimates  form  the  basis  ol  his  estimated  failure  distribution  F(t) 
Using  this  in  place  of  F(t^)  In  (3)  and  minimizing  with  respect  to  t^ 

A Jfc 

produces  t^  from  which  one  can  compute 

G(e)  - U(t£)  - U*  , (repeated  here  from  (8))  . 

This  procedure  was  repeated  1000  times  for  each  n and  then  repli- 
cated. Examples  of  distribution-*  of  G(n)  are  shown  In  graphs  10  and 
11  In  the  Appendix.  The  results  were  averaged  obtaining  point  estimates 


* It  was  noted  that  the  distribution  of  G(n)  is  well  fit  by  a Weibull 
distribution.  The  variance  of  the  distribution  Increases  with  the 
parameter  A . 


of  G(n)  « A smooth  curve  vss  passed  thru  these  points  using  the 
function 

G<n)  ■ exp[  s + b £n(n)]  . 

The  result  is  shown  plotted  in  graph  5.  Finally  forward  differences 
in  G(n)  were  taken,  A G(n)  . These  are  shown  in  graph  6 for 
A ■ 100  and  in  graph  7 for  A * 200  . As  an  illustration  of  the 
utility  of  these  graphs  in  finding  n * , wa  take  the  following 

example) 

tasting  cost/tube  • $300,000 
cost  of  tuba  » $25,000 

total  expected  utilization  of  the  system  • 1 million  rnds. 

Thus,  B • 300,000/25,000  - 12  ind  106  B/u  • 12  . 

From  graph  5,  we  see  that  10^  A G(n)  > 12  for  n » 18  and 

106A  G(n)  < 12  for  n - 19 

Therefore  a*  • 18  • 
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SUMMARY 


Two  problems  are  addressed  In  this  paper:  determination  of  the 
optional  tube  safe-life  and  determination  of  the  optimum  size  of  sample 
to  be  life  tested.  These  problems  are  related  in  that  the  method  by 
which  the  decision-maker  would  choose  a "safe"  life  after  tests  are 
completed  affects  the  penalty  associated  with  imperfect  information. 

This  had  a bearing  on  the  number  of  samples  taken.  Stated  succinctly, 
the  way  in  which  the  information  is  used  affects  the  amount  of  effort 
productively  spent  in  its  acquisition.  In  an  obvious  way.  the  informa- 
tion concerning  failure—  once  acquired— should  be  used  with  all  of  the 
pertinent  cost  factors  to  choose  an  optimum  safe-life.  In  this  note 
we  have  shown  how  this  can  be  dona.  We  have  also  shown  that  three 
principle  parameters  srs  involved  in  these  decisions:  A,  B,  and  u . 
further,  since  the  decision-making  models  Involve  futore  (and  uncertain) 
variables,  vs  have  described  the  sensitivity  of  the  decision  to  the 
parameters. 

Finally  one  might  conjecture  that  the  decisions  themselves  affect 
the  future  values  of  some  of  the  variables.  For  example,  a decision  to 
lower  the  safe-life  raises  the  total  inventory  requirement.  In  an 
attempt  to  meet  the  added  requirement  ueege  may  decrease  due  to  rationing 
of  tubes.  -In  some  caeca  this  reciprocal  relationship  between  e decision 
and  its  data  base  may  not  be  negligible. 


COMPUTER  PROGRAMS 


The  computer  program  implementation  of  the  analyai8  is  simple 
and  direct*  Basically,  the  program  finds  the  expectation  of  G(n) , (J(n), 
for  values  of  n from  a set  read. in  as  parameter  data.  Other  values 
treated  as  parameters  are  El,  E2,  E3  and  A.  El,  E2  and  E3  are, 
respectively,  values  of  t in  ?(t)  for  which  F(t)  - .16,  .30,  and 
.84  . Thus,  for  ?(t)  a log-normal  distribution  with  parameters 
u and  a , u • in  E2  and  a • 1/2  In  (E3/E1)  . 

The  program  was  written  in  the  FORTRAN-! I language  for  implemen- 
tation on  a FORTRAN-17  compiler.  The  program  usee  several  standard 
library  programs  developed  by  this  office.  Since  these  programs  are 
not  part  of  a basic  FORTRAN  system,  their  source  listings  will  appear 
in  this  note  for  completeness.  The  programs  used  are 
(1)  SUBROUTINE  MIN  10  (X0,  70,  Y,  DX,  EPX,  EPY) 

Y is  a function  name  which  is  given  in  the  EXTERNAL 
statement  of  the  calling  program.  The  program  minimizes 
Y(X),  a user-writ tan  function,  by  the  method  of  interval- 
halving. A starting  value  X0  , hopefully  near  the  optimum, 
ia  given.  The  steps  in  X of  DX  are  used  until  the 
aearch  is  terminated  whan  changes  in  both  X and  Y 
ltarataa  era  less  than  the  given  epelloo  values  EPX  and 
EPY  , respectively.  The  optimum  point  Is  passed  back  to 
the  calling  program  thru  X0  and  Y0. 
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(2)  SUBROUTINE  OPEN  (IHIST,  XMIN,  INTER,  DX,  Dl,  D2,  D3,  XMAX) 
SUBROUTINE  DIST  (IHIST,  XMIN,  INTER,  DX,  Dl,  D2,  D3,  X) 
SUBROUTINE  CLOSE  (IHIST,  XMIN,  INTER,  DX,  Dl,  D2,  D3,  N) 

Thii  sec  of  programs  is  used  to  construct  a histogram 
for  the  FORTRAN  random  variable  X . IHIST  Is  the 
name  of  a one-dimensional  array  used  for  storing  the 
frequency-of-occurrenca  of  X . XMIN  Is  the  lover 
bound  of  the  histogram,  INTER  is  the  number  of  intervals 
in  the  histogram.  DX  is  the  interval  width.  Dl,  D2  and 
D3  ere  dummy  variables  used  by  the  subroutines*  Note  that 
DX,  Dl,  D2  and  D3  ere  computed  and  used  only  by  the  hi • co- 
gram  package.  XMAX  is  the  upper  limit  of  the  histogram. 

X is  the  random  variable  to  be  DlST-rlbutad  and  N is 
the  number  of  random  variables  that  were  distributed  for 
final  summary  computations  in  CLOSE. 

(3)  FUNCTION  SNORM  (Z) 

The  program  computes  the  Integral  from  — to  Z of 
the  standard-normal  probability  function.  The  error, 

/c(x)  / <,  1.5  x 10  . 

(4)  RANDOM  (N) 

This  function  subprogram  was  written  in  HA?  language  for 
the  IBM-7044.  Each  time  it  appears  and  Is  called,  a ran- 
dom number  is  produced.  That  number  is  the  sum  of  N 
uniform,  random  deviates  from  (0,1)  . It  is  used  in 
this  study  to  generate  normal  deviates  using  the  Central 
Limit  thaortm  technique  as  explained  in  OR  Technical  Note 


66-1  by  Stuart  Olson. 
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The  main  program  begins  by  rsading  a aet  of  constants,  Identi- 
fication and  paraaatara.  The  FORI RAH  symbols  sad  thair  meaning  are 
given, 

Symbol  Meaning 

TITLE  An  array  used  to  store  80  characters  of  run 

Identification, 

MS  Number  of  Iterations  to  be  used  in  Monte  Carlo 

simulation  (1000). 

NPAR  Humber  of  parameter  cards  in  the  set. 

IDL  A number  changed  for  each  run  so  that  the  random 

number  sequence  does  not  start  with  the  seme  number. 
HVHT  Humber  of  intervals  In  the  G(n)  histogram  XVN. 

HINT  Number  of  intervale  In  the  e£  histogram  ITS. 

HUNT  Number  of  intervals  in  the  U(t£)  histogram  XUS. 

TS6  An  estimate  of  (1100). 

OT  An  Increment  used  In  the  search  for  optimum  t*  (25). 

Is 

ElgEU  Mixlsuai  fraction*!  error  Uvtla  for  U(t  ) 

is 

mlnimixetlao  used  by  HXNlD  (.002  and  .001). 

VL,VH  Lowsr  end  upper  limits,  respectively,  for  G(n) 

histogram  (0.0  and  0.01). 

TL,TH  Lower  end  upper  limits,  respectively,  for  t* 

la 

histogram  (0.Q  and  10000.). 

UL.UH  Lovsr  end  upper  limits,  respectively,  for  U(l*) 

la 

histogram  (0.U01  end  0.005). 
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(The  following  are  repeated  NPAR  times.) 

AV(I)  i-th  value  of  A in  the  parameter  set 

B1V(I) ,E2V(I) ,E3V(I)  i-th  values  of  El,  E2  and  E3 

(2600,  4050  and  6300,  for  all  1 ) 

NV(I,J)  For  the  i-th  set  of  perameters,  the  values 

of  n to  be  used 

Using  the  parameters  given,  the  program  computes  y and  o . 
These  nre  used  to  find  t*  by  the  MIN  ID  search.  A loop  counter  is 

U 

initialized  to  sweep  thru  ail  the  n's  given  on  the  i-th  parameter 
card.  The  histogram  files  are  OPEN-ed  and  the  Monte  Carlo  loop  is 
initialized.  Within  the  loop,  y and  a are  obtained  from  equation 
(14) < When  these  are  computed,  the  value  of  tL  must  be  found  such 
that  U(tL)  is  a minimum.  This  is  accomplished  by  a special  search 
subroutine  called  MINFH  which  searches  on  the  variable 
z • (in  tT  - y )/a  . When  t£  is  found,  GN  » U(t*)  - U(t*)  is 
computed  and  distributed  in  the  IVN  histogram.  Other  interesting 
quantities  are  also  distributed  in  other  histograms,  After  NS  repeti- 
tions, the  current  operating  paramatars  are  printed  and  the  hlatogram 
filet  are  ClOSE-d  end  printed.  Then  the  n-loop  is  repeated,  followed 
by  the  NPAR  loop. 

The  progrem  source  listings,  sample  input  and  output  are  given. 
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r/i  O O 


siesYs 

SJOB 

STIHE 

SPAGfc 


SOURCE  LISTINlOS 

AWCOR-SQ 

NC 

300 


IBJO0  UCC  HAP i DECK 
IBFTC  GUNSS 


JUNE  5,196? 

GUN  TUBE  SAMPLE  SIZE  STUDY  SCHLENK6R  - OLSON 

extehnal  u 

DIMENSION  AV[2Q]#ElVt2Q]»62Vt20)»E3Vl2Q1»NVt20,10I,TITLEtl6]# 

1 I VN  (10021  , ITS  l5021,,IUSI2fl2  1,12  [102  J 
COMMUN  A,EHU,SIQ,EMUH,SI»iH,TLO,UHS,Z,ZSTAR 
READ  PARAMETERS  FOR  SIMULATIONS 

READ  i»title,ns«npar,iul,nvnt,ntnt,nunt»tso,dt,et,eu,vl,vm,tl, 

1 TH.UL.UM.  [AV(1 1 ,E1V(  1 1 ,62V(1 3 , E3V  (II,  (NVU,  J]  , Jsl,10I  , I*1,NPAR] 
1 FORMAT  C16A5/I5. 12, 13, 3 1 5 , 5X4F1 0 . 0 /6F 10 . 0 / 1 4F 1 0 . 0.1012)3 
C IDLE  ranuoh  NUMBER  GENERATOR  TO  start  NEW  SEQUENCE 
DO  9 1=1, ICL 
9 HANDsRANDOM(l) 

C SET  LIMITS  OF  2 HISTOGRAM 


ZL  = -3.5 


ZH=-2.5 

C INITIALIZE  PARAMETER-SET  COUNTER 

DO  2 I P A R = 1 , NPAR 

C GET  WORKING  VALUES  OF  PARAMETERS  FROM  HEAD-IN  ARRAY 

AsAVCIPAR] 
tl =ElV ( I PAH  1 
E2*E2V[ I PAR ] 


E3=E3VI IRAK] 

C COMPUTE  MU  AND  SIGMA  FOR  UITLI  FUNCTION 
fcMU=AL0Gl&2 J 
S I G* . 5* AlOG  t E 3/61 5 

C COMPUTE  SIGMA  WITH  CORRECTION  FACTOR  FOR  VARIANCE  IN 

C CENTRAL  LIMIT  SAMPLING  SCHEME  to  be  used  below 

blliCs.  70/10678-SIG 

C INITIALIZE  tl  TO  86  USED  IN  SEARCH  FOR  MIN-U  ' 

TlSaTSO 

C COMPUTE  EPSILON  VALUeS  TO  BE  USED  IN  GENERAL  ONE-DIMENSIONAL 

C SEARCH  FOR  MJN-U 

EPT=6T.TUS 
EPU=EU*U(TLS1 

c call  ONE-OJMENSJONAI  SEARCH  TQ  MINIMIZE  U(Tll  FUNCTION  AND 

C STuRs  RESULT  IN  TLS  AND  USTAR 

CALL  MIN1D(TLS,USTAR,U.DT,£PT.6PU) 

ZSTAH*( AlOGITLSI -EMUi/SIG 

C PRINT  RESULTS  OF  SEARCH  PLUS  A TABLE  OF  U(Tl) 

PRINT  11 » TLS , UbT  AR 

11  FOHHaT  (JlHlTABLt  OF  UITLI.  MINIMUM  AT  TL»F8 ♦ 2, 3X7HU-STAH*F U . 6/ 
1 1W013X2MTL10X5HUITU  I 
TLT*Q  . 
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DO  12  1*1*50 

TLT*TLT  + 10Q , • . 7 7 , . 

UTL*UCTLT]  ) 

12  PRINT  13*  TLT » UTL 

13  FORMAT  tlH  F15,2*F15,8) 

C INITIALIZE  COUNTER  FOR  N-VALUES 

DO  2 NN*1 * 10 
N = NV ( I PAR*  NN  3 

C PERFORM  SIMULATION  USING  N-VALUE  GIVEN 

IF  ( N-l 3 2 * 2 * 3 

3 FN*N 

C OPEN  HISTOGRAM  DATA-FILES 

CALL  OPEN!  I Z , ZL*  100 , 241,  Z42 . Z4'3*  ZHJ 
CALL  OPEN!  1VN.VL, NVNT,ZU,Z12* 213.  VH) 

CALL  OPEN! ITS* TL , NTNT , Z21 , 222  * Z23 » TH 1 
CALL  OPEN  I IUS»UL»NUNT*Z31,Z32*Z33*UH) 

C INITIALIZE  SIMULATION  COUNTER 

DO  4 I si # NS 

C COMPUTE  MU-HAT  AND  SIGMA-HAT  VALUES 

kMUHsfl. 

SIGHaQ. 

DO  5 J=1*N 

XsEMU*SIGC*CRANOOM(243 -12. 3 
EMUHagMUH+X 

5 S I GHsS I GH*X*X 
tMUHsEMUH/FN 

SIGHsSQRT  USIGH-FN*fcMUH#6MUH)/tFN-l.  1 J 
C USING  EMUH  AND  SIGH,  MINIMIZE  UITL3  FUNCTION 

CALL  MJNFH 

C USING  EMU  AND  SIG  VALUES  WITH  THE  OPTIMUM  TL  VALUE.  TLO,  JUST  ) 

C OBTAINED,  COMPUTE  UITL03 

UTS=U I TLO 1 
GN=UTS-USTAP 

C 0ISTHI8UTE  COMPUTED  VALUES  IN  THEIR  RESPECTIVE  HISTOGRAMS 
CALL  DISTl iVN.VL,NVNT,Zll.Z12.Z13,GN3 
CALL  DISTI  IZ.ZL * 10  0.241, Z42,r43*Z 3 
CALL  DlSTt ITS. TL, NTNT, Z2l.Z22.Z23, TLO) 

4 CALL  0ISTC IUS.UL,NUNT,Z31,Z32.Z33*UTS1 

C END  OF  A SIMULATION-PRINT  RESULTS 

PRINT  6*  T1 TLE* A . El , 62 . 63 , EMU , S I G. N 

6 FORMAT  t 1H116A5/1M014X1MA13X2HS113X2WE213X2H63/1H  4Fl5 , 2/1H013X 
1 2HMU10X3HSIGMA14X1MN/1H  2F15.4.U5) 

PRINT  14 

14  FORMAT  I16H1GIN)  STATIST  I CS  J 

C CLOSE  HISTOGRAM  FILES  AND  PRINT  RESULTS 
CALL  CLOSE  tlVN,VL.NVNT.m,Zl2.Zl3.  NS  J 
PRINT  7 

7 FORMAT  (19H1TL-STAR  STATISTICS! 

CALL  CLOSE  HTS.TL.  NTNT,  Z21.Z22.Z23,  NS  J 
PRINT  8 
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6 FORMAT  {18H1U-STAR  STATISTICS] 

CALL  CLOSE! IUS,UL»NUNT»Z31>Z32.Z33,N$1 
PRINT  10 

10  FORMAT  C13H1Z-STATISTICS) 

CALL  CLOSE  UZ.ZL,  100. Z41.Z42.Z43, NS J 
2 CONTINUE 
CALL  EXIT 


SI8FTC  UTL 

FUNCTION  U £ TL  3 
COMMON  A * EMU  > S IG 

C SNQRM  computes  THE  STANDARD-NORMAL  INTEGRAL  FOR  ANY  Z VALUE.  IT  IS 
C ACCURATE  TO  1.5  PARTS  IN  THE  7-TH  OECIMAL 
Us(A*SnORM([ALOG{TL1-EMU]/SI61*1.]/TL 
RETURN 
END 


; 
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SI8FTC  M I NFH 

C THIS  ROUTINE  MINIMIZES  UITL3  BY  SEARCHING  ON  THE  VARIABLE  Z 
SUBROUTINE  M I NFH 

COMMON  A.EMU,SIG,EMUH,SIQM,T2,U2»ZiZSTAR 
Z=ZSTAR 

T = EXP  [ S I GH-Z+EMUH 3 
U*[A*SNORMIZJ*l.I/T 
DZ«.l 

8 Z2*Z*DZ 

T2*EXP  (SIGH*Z2*EMUHJ 
F sSNORM ( Z2  3 
U2«tA*F*l. J/T2 
IF  [ABS  [T-T23-1.034.4.5 

4 IF  (APS  CU-U2J-U*!. 6-319. 9.5 

5 IF  (U2-U36.6.7 
7 0Z*-.5*DZ 

6 USU2 
1-12 
T = T2 

GO  TO  8 

9 RETURN 
END 


SIBFTC  MIN1D 

C MJNIU  IS  A SIMPLE  QNE-QIMENS I ONAL  FUNCTION MINIMI 2ER 
SUBROUTINE  M I NlD l XS » YS> FUNC # DXS, EFX, EPY ) 

x=xs 

Y = FUNC I X ] 

UX=DXS 
105  X2sX+DX 

Y28FUNCCX2I 

IF  (ABS  CDX!-EPX)10i,101,102 

101  IF  (ABS  (Y-Y2)-EPYJ100#100.102 

102  IF  I Y2- Y 1 1 Q3  # 104  » 104 
104  0X=- . 5*DX 

103  YsY2 
' X = X2 

GO  TO  105 
100  XS=X2 
YS  = Y2 
RETURN 
END 


SIbFTC  OPEN 

SUBROUTINE  OPEN  UP, XM IN, JNT , DX » DS »DSSQ, XMAX ) 

DIMENSION  I P 1 2502 1 

DSsO. 

DSSQsO . 

LIMsJNT*2 
DO  1 Isl,LlM 

1 ipm  = o 

DXs  [XMAX-XM1NI/FL.QAT  UNT] 

RETURN 

END 
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S 1 BFTC  D I ST 

SUBROUTINE  DISTtIP,XMlN#JNT,DX,DS,QSSQ,.X] 

DIMENSION  I P 1 2502 3 

DS=DS+X 

USSO=DSSQ^X*X 

Y = X - XM  1 N 

KsY/OX^.5 

IF  C K 3 1#2»3 

3 IF [K-JNT ] 6>  4<5 

6 IF  tY-FUOAT  I K ] *DX 1 2*4*2 
2 KsK*1- 

4 IPtK]sJP[K]*l 
HETURN 

1 K= JNT+1 
GO  TO  4 

5 KsjNT+2 
GO  TO  4 


SlbFTC  CLOSE 

SUBROUTINE  CLOSE  1 1 Pi XMIN, I NTER . DX , SS, SSSS . N 1 
DIMENSION  I P l 25 02 3 
FN  = N 

DS=S$/FN 

DSSO=SURT  (ABS  U SSSS-FN*DS*DS J / C FN-1 . J 1) 

SX  = 0 . 

SIUX50  . 

V=XMlN-DX/2. 

UO  1 1=1# INTER 
' V = V + DX 

DEN  = FLQAT  [IPim*V 
SX=SX*UEN 
1 j>IUX  = S IGX*UEN*V 
SX=SX/FN 

S I liXsSUWT  t A8S  [SIi»X/FN-SX*SX-UX*DX/12»ll 
PRINT  2,N,0S,DSSU,SX,SIGX 

FORMAT  [1H010X26HHISTOGRAM  GENERATOR  OUTPUT  /1H  9X6HSAMPLE6X 

1 9HSAMP  MEAN3X12HSAMP  STD  DEV6X9HHIST  .NEAN3X 12HH I ST  STU  DEV/1H 

2 I15,1P4c1S.6/9H0INTERVAL4X11HUPPER  L 1 N I T5X10HCUSS  MARK3X4HFRE08X 

3 7HDtNSI TYbXIGHCUMULATl VE3 
L I nE=»40 

UL  s X M I N 

CUMsFLOAT  l IP C INTER*1]]/FN 
PRINT  H, JL# I P C INTER-1] *CUM,CUM 
8 FORMAT  tlH  7XlM01Ptl5.4,15Xl7.0P2Fl5.SJ 
UMsUL  + DX /2 « 

UL=UL-DX 

DO  3 1=1, INTER 

IF  t IP( I ] 110,11# 10 

10  UEN=FLUAT  l IP( I Il/FN 
CUM=CUM*UfcN 

LlNE  = UNs*l 
IF  t L I NE ) S » 5 , 6 
a IF  t !NT6W-l-05i5*5»l5 
la  * = 0 

DO  12  Js I, INTER 
IF  1 1 P t J ) )1?, 12,13 
13  Ksk*1 
12  CONTINUE 

IF  CK-0515.5,14 
1**  PRINT  7 

7 FORMAT  ( lHl8HlNTfcRvAL4XHMUPP£R  L I H I T5X1 0«CL ASS  HARK3X4HFRfcQ8X 
1 7HDtNslTY5XiQHCUMULATlve> 

Li NEs-55  . 

a PRINT  4, 1 .UL.CM, I P ( I J , DEN, CUM 
h FORMAT  (1M  I8,1P2E15.4,I7,0P2F15.8) 

11  JL=Ul*DX 
o Cm sCM* D X 

K= |NTER*2 


I =K-1 

DEN=FLQAT  (IPIKll/FN 
CUM  = CUM  + L)£N 

PRINT  9, I*  I P t K J , DEN, CUM 
9 FOHMAT  C1H  I 8 » 30 X I 7 , 2F 15 . 8 ] 
KETUPN 
END 


SIbFTC  SNOkM 

FUNCTION  SNORMIZ) 

IF  (APS  IZ1-7. 11.1,2 

1 SNORM=. 5-SIGN  [ .5,Z3*ERRFX[ABS  (ZJ*. 707106781 
RETURN 

2 SNQR Ms. 5-SIGN  1.5,1) 

RETURN 

END 

SItiFTC  ERRFX 

FUNCTION  ERRF  X [ X ] 

Tsl./tl.*.3275911*Xl 

SUM=1. 0614Q543*T 

SUM= [SUM-1 • ^5  31520  3 3 *T 

SUHs(SUM*1.421413741*T 

SUM=CSUM-.284496736]*T 

iUM=(SUM*.254829592!*T 

tRRF X*i .-SUM*EXP  t-X*X] 

kETURN 

fcND 


SIBHAP  RANDOM 
ENTRY 

RANDOM  SAVE 
TZE 
TPL 
CHS 
STO 

RETURN 
NEWUO  CLA 

RETURN 
START  CLA 
STO 
CLA* 

gener  pax 

LOO 

ZAC 

MPY 

LLS 

ZAC 

LLS 

STO 

OkA 

PAD 

PAD 

STO 

TIX 

RcTURN 

sum  ess 

UO  OCT 

MASK  OCT 

X DEC 

END 


RANDOM 

1.2.4 
NEWUO 
START 

UO 

RANDOM 

UO 

RANDOM 
= 0 
SUM 

3.4 
.1 
UU 

X 

8 

27 

UO 

MASK 

MASK 

SUM 

SUM 

GENER* 1 #1.1 
RANDOM 
1 

00042-3071643 
200000000000 
1 15d  7 
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SAMPLE  INPUT 


01234567b9Q1234567d90l2345678901234567890l234567B90l2345676901234567890l23456789  COL.  SC  AlE 


UON  T USE 

Sample 

SIZE  STUDY 

ANC-OR 

1000  i 

11  10«0 

500  200 

1100. 

0. 

.01 

.0 

10000 

200  . 

2600  . 

4050. 

6300  . 

<00. 

2600  . 

4050. 

6300  . 

100. 

2600  . 

4050  . 

63u0  . 

SCHLENKbR  * OLSON  JUNE. 20. 1967 
25.  .00?  .001 

0.001  0.005 

81012141616202224') 
81012141618202224  ? WrA**3 
81012141618202224^ 


TABLE  OF  U (TL) • MINIMUM  AT  TL*  989.06  U-STAR*  0.00115709 


TL 

U(TL ) 

100.00 

0.01000000 

200.00 

0.00500000 

300.00 

0.00333334 

400.00 

0.00250004 

500.00 

0.00200046 

600.00 

0.00166933 

700.00 

0.00143898 

800,00  . 

0.00128092 

900.00 

0.00118630 

1000.00 

0.00115735 

1100.00 

0.00120230 

1200.00 

0.00133180 

1300.00 

0.00155625 

1400.00 

0.00188393 

1500.00 

0.00231984 

1600.00 

0.00286522 

1700.00 

0.00351750 

1800.00 

0.00427067 

1900.00 

0 .0051 1584 

2000.00 

0.00604189 

2100.00 

0.00703621 

2200.00 

0.00808532 

2300.00 

0.00917547 

2400.00 

0.01029313 

2500.00 

0.01142533 

2600.00 

0.01256002 

2700.00 

0.01368620 

2800.00 

0.01479408 

2900.00 

0.01587513 

3000.00 

0.01692206 

3100.00 

0.01792882 

3200.00 

0.01889055 

3300.00 

0.01980343 

3400.00 

0.02066467 

3500.00 

0.02147235 

3600.00 

0,02222534 

3700.00 

0.02292321 

3800.00 

0.02356613 

3900.00 

0.02415473 

4000.00 

0.02469011 

4100.00 

0.02517368 

4200.00 

0.02560712 

4300.00 

0.02599234 

4400.00 

0,02633139 

4500.00 

0.02662646 

4600.00 

0.0268797? 

4700.00 

0.02709362 

4800.00 

0.02727029 

4900.00 

0.02741205 

5000.00 

0.02752115 

GUN  TUBE  SAMPLE  SIZE  STUDY  AWC-OR  SCHLENKER  - OLSON  JUNEt20,1967 


A 

200.00 


El 

2600.00 


E2 

4050.00 


MU 

8.3065 


SIGMA 

0.4425 


N 

8 


E3 

6300.00 
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G { N ) STATISTICS 


HISTOGRAM  GENERATOR  OUTPUT 


SAMPLE  SAMP  MEAN  SAMP  STD 

1000  8.122838E-04  1.703063E 

DEV 

-03 

HIST  MEAN 
7.437439E-04 

HIST  STD  DEV 
1.3602276-03 

f AL 
0 

UPPER  LIMIT 

0. 

CLASS  MARK  FREQ 

0 

OENSITY 

0.00000000 

CUMULATIVE 

0.00000000 

l 

1 .00006-05 

5.00006-06 

73 

0.07300000 

0.07300000 

2 

2 .0000E-05 

1.50006-05 

60 

0.06000000 

0.13300000 

3 

3.0000E-05 

2.5000E-05 

38 

0.03800000 

0.17100000 

A 

4.00006-05 

3.50006-05 

27 

0.02700000 

0.19800000 

5 

5.0000E-05 

4.50006-05  ’ 

31 

0.03100000 

0.22900000 

6 

6.Q000E-05 

5.50006-05 

19 

0.01900000 

0.24800000 

7 

7.00006-05 

6. 5000E-05 

18 

0.01800000 

0.26600000 

8 

8.0000E-05 

7.50006-05 

21 

0.02100000 

0.28699999 

9 

9.00006-05 

8.50006-05 

24 

0.02400000 

0.31099999 

10 

1.00006-04 

9.50006-05 

21 

0.02100000 

0.33199999 

ll 

1.1000E-04 

1 • 0500E-04 

14 

0.01400000 

0.34599999 

12 

1.2000E-04 

1 . 1500E-04 

20 

0.02000000 

0.36599999 

13 

1.3000E-04 

1.25006-04 

9 

0.00900000 

0.37499999 

14 

1.4000E-04 

1.35006-04 

13 

0.01300000 

0.38799998 

IS 

1.5000E-04 

1.45006-04 

17 

0.01700000 

0.40499998 

16 

1.6000E-04 

1.5500E-04 

15 

0.01500000 

0.41999998 

17 

1.70006-04 

1.6500E-04 

12 

0.01200000 

0.43199997 

18 

l .80006-04 

1.75006-04 

14 

0.01400000 

0.44599997 

19 

1.9000E-04 

1.85006-04 

9 

0.00900000 

0.45499997 

20 

2 .00006-04 

1.9500E-04 

9 

0.00900000 

0.46399997 

21 

2.10006-04 

2.0500E-04 

16 

0.01600000 

0.47999997 

22 

2. 20006-04 

2 . 1500E-04 

9 

0.00900000 

0.48899997 

23 

2.30006-04 

2.25006-04 

10 

0.01000000 

0.49899997 

24 

2 .40006-04 

2.35006-04 

ll 

0.01100000 

0.50999997 

25 

2.5000E-04 

2.4500E-04 

13 

0.01300000 

0.52299996 

26 

2.60006-04 

2,55006-04 

9 

0.00900000 

0.53199996 

27 

2.7000E-04 

2.6500E-04 

6 

0.00600000 

0.53799996 

28 

2.80006-04 

2.75006-04 

2 

0.00200000 

0.53999995 

29 

2.90006-04 

2.85006-04 

9 

0.00900000 

0.54899995 

30 

3.00006-04 

2.95006-04 

7 

0.00700000 

0.55599995 

31 

3.10006-04 

3.05006-04 

11 

0.01100000 

0.56699995 

32 

3.20006-04 

3.15006-04 

3 

0.00300000 

0.56999995 

33 

3.30006-04 

3.25006-04 

6 

0.00600000 

0.57599995 

34 

3. 40006-04 

3.35006-04 

3 

0.00300000 

0.578999V4 

35 

3.50006-04 

3.45006-04 

7 

0.00700000 

0.58599994 

36 

3.60006-04 

3.55006-04 

7 

0.00700000 

0.59299994 

37 

3.70006-04 

3.65006-04 

6 

0.00600000 

0.59899994 

38 

3. 80006-04 

3.75006-04 

7 

0.00700000 

0.60599994 

39 

3.90006-04 

3.85006-04 

11 

O.OilOOOOO 

0.61699994 

40 

4.00006-04 

3.9500E-04 

5 

0.00500000 

0.62199993 

37 


NTERVAt 

UPPER  LIMIT 

CLASS  MARK 

41 

4.1000E-04 

4. 05006-04 

42 

4.2000E-04 

4.1500E-04 

43 

4. 3000E-04 

4.2500E-04 

44 

4.4000E-04 

4.3500E-04 

45 

4.5000E-04 

4.4500E-04 

46 

4.6000E-04 

4.55006-04 

47 

4.7000E-04 

4.6500E-04 

48 

4.8000E-04 

4. 75006-04 

49 

4.9000E-04 

4.85006-04 

50 

5.0000E-04 

4.9500E-04 

51 

5.1000F-04 

5.0500E-04 

52 

5.2000E-04 

5.1500E-04 

53 

5. 3000E-04 

5. 2500E-04 

54 

5.4000E-04 

5.3500E-04 

55 

5.5000E-04 

5.4500E-04 

56 

5.6000E-04 

5.5500E-04 

57 

5.7000E-04 

5.6500E-04 

58 

5.8000E-04 

5.7500E-04 

60 

6.0000E-04 

5.9500E-04 

61 

6. 1000E-04 

6. 0500E-04 

62 

6.2000E-04 

6.1500E-04 

63 

6.3000E-04 

6.2500E-04 

64 

6.4000E-04 

6* 3500E-04 

65 

6.5000E-04 

6.4500E-04 

66 

6.6000E-04 

6. 5500E-04 

68 

6.8000E-04 

6.7500E-04 

69 

6.90006-04 

6.8500E-04 

70 

7.0000E-04 

6.9500E-04 

71 

7.1000E-04 

7«0500E-04 

72 

7.2000E-04 

7.1500E-04 

73 

7.3000E-04 

7.2500E-04 

74 

7.4000E-04 

7.3500E-04 

75 

7.5000E-04 

7.4500E-04 

76 

7.6000E-04 

7.5500E-04 

77 

7.7000E-04 

7.6500E-04 

78 

7.8000E-04 

7.7500E-04 

79 

7.9000E-04 

7.8500E-04 

80 

8.0000E-04 

7.9500E-04 

81 

6.10006-04 

8.0500E-04 

82 

8.2000E-04 

8.15006-04 

83 

8. 300QE-04 

8.2500E-04 

84 

8.4000E-04 

8. 3500E-04 

85 

8.5000E-04 

8. 4500E-04 

86 

8.6000E-04 

8.5500E-04 

88 

8.8000E-04 

8.7500E-04 

89 

8.9000E-04 

8.8500E-04 

90 

9.0000E-04 

8.9500E-04 

91 

9. 1000E-04 

9.0500E-04 

92 

9. 2000E-04 

9. 1500E-04 

93 

9. 3000E-04 

9. 2500E-04 

94 

9.40006-04 

9. 3500E-04 

95 

9.5000E-04 

9.4500E-04 

96 

9.6000E-04 

9.5500E-04 

98 

9.8000E-04 

9. 7500E-04 

100 

1.0000E-0J 

9. 9500E-04 

101 

1.0100E-03 

1.0050E-03 

FREQ 

DENSITY 

CUMULATIVE 

4 

0.00400000 

0.62599993 

6 

0.00600000 

0.63199992 

5 

0.00500000 

0.63699992 

4 - 

0.00400000 

0.64099991 

2 

0.00200000 

0.64299991 

4 

0.00400000 

0.64699990 

5 

0.00500000 

0.65199990 

3 

0.00300000 

0.65499990 

4 

0.00400000 

0.65899989 

4 

0.00400000 

0.66299988 

5 

0.00500000 

0.66799988 

4 

0.00400000 

0.67199987 

2 

0.00200000 

0.67399987 

3 

0.00300000 

0.67699987 

7 

0.00700000 

0.68399987 

4 

0.00400000 

0.68799986 

1 

0.00100000 

0.68899985 

1 

0.00100000 

0.68999985 

3 

0.00300000 

0.69299985 

2 

0.00200000 

0.69499984 

2 

0.00200000 

0.69699984 

5 

0.00500000 

0.70199984 

6 

0.00600000 

0.70799983 

1 

0.00100000 

0.70899983 

5 

0.00500000 

0.71399982 

3 

0.00300000 

0.71699982 

2 

0.00200000 

0.71899982 

5 

0.00500000 

0.72399981 

5 

0.00500000 

0.72899981 

3 

0.00300000 

0.73199981 

1 

0.00100000 

0.73299980 

2 

0.00200000 

0.73499980 

4 

0.00400000 

0.73899979 

4 

0.00400000 

0.74299978 

3 

0.00300000 

0.74599978 

1 

0.00100000 

0.74699978 

2 

0.00200000 

0.74899977 

2 

0.00200000 

0.75099977 

1 

0.00100000 

0.75199977 

2 

0.00200000 

0.75399976 

3 

0.00300000 

0.75699976 

3 

0.00300000 

0.75999976 

2 

0.00200000 

0,76199976 

1 

0.00100000 

0.76299975 

2 

0.00200000 

0.76499975 

l 

0.00100000 

0.76599974 

2 

0.00200000 

0.76799974 

4 

0.00400000 

0.77199973 

4 

0.00400000 

0.77599972 

1 

0.00100000 

0.77699972 

2 

0.00200000 

0.77999972 

l 

0.00100000 

0.77999971 

2 

0.00200000 

0.78199971 

2 

0.00200090 

0.78399970 

2 

0.00200000 

0.78599970 

3 

0.00300000 

0.78899970 

38 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

FREQ 

DENSITY 

CUMULATIVE 

102 

1.02006-03 

1.0150E-03 

3 

0.00300000 

0.79199970 

103 

1.03006-03 

1.0250E-03 

1 

0.00100000 

C. 79299969 

104 

1.04006-03 

1.03506-03 

3 

0.00300000 

0.79599969 

105 

t .05006-03 

1.04506-03 

1 

0.00100000 

0.79699969 

107 

1.07006-03 

1.06506-03 

2 

0.00200000 

0.79899968 

108 

1.08006-03 

1.07506-03 

l 

0.00100000 

0.79999968 

109 

1.09006-03 

1.0850E-03 

2 

0.00200000 

0.80199967 

111 

1.11006-03 

1 • 1050E-03 

1 

0.00100000 

0.80299967 

112 

t. 12006-03 

l . 1 1 50E-03 

3 

0.00300000 

0.80599967 

113 

1.13006-03 

1.12506-03 

. 3 

0.00300000 

0.80899967 

115 

1.15006-03 

1.14506-03 

l 

0.00100000 

0.80999966 

116 

1.16006-03 

1.15506-03 

3 

0.00300000 

0.81299966 

117 

1.17006-03 

1.16506-03 

2 

0.00200000 

0.81499965 

118 

1.18006-03 

1 • 1750E-03 

1 

0.00100000 

0.81599965 

120 

1.20006-03 

1.19506-03 

l 

0.00100000 

0.81699964 

121 

1.21006-03 

1.20506-03 

5 

0.00500000 

0.82199964 

122 

1.22006-03 

1.21506-03 

3 

0.00300000 

0.82499964 

124 

1.24006-03 

1.23506-03 

2 

0.00200000 

0.82699963 

126 

1.26006-03 

1.25506-03 

1 

0.00100000 

0. 82799963 

130 

1.30006-03 

1.29506-03 

1 

0.00100000 

0.82899962 

132 

1.32006-03 

1.31506-03 

1 

0.00100000 

0.82999962 

133 

1.33006-03 

1.32506-03 

1 

0.00100000 

0.83099961 

134 

1.34006-03 

1.33506-03 

l 

0.00100000 

0.83199961 

135 

1.3500E-03 

1.34506-03 

2 

0.00200000 

0.83399960 

136 

1.36006-03 

1.35506-03 

1 

0.00100000 

0.83499960 

137 

1.37006-03 

1.3650E-03 

l 

0.00100000 

0.83599959 

138 

1.38006-03 

1.37506-03 

1 

0.00100000 

0.83699959 

139 

1.39006-03 

1.38506-03 

l 

0.00100000 

0.83799958 

140 

1.4000E-03 

1.39506-03 

1 

0.00100000 

0.83899958 

141 

1.41006-03 

1.40506-03 

l 

0.00100000 

0.83999957 

142 

1.42006-03 

1.41506-03 

1 

0.00100000 

0.84099957 

144 

1.44006-03 

1.43506-03 

1 

0.00100000 

0.84199956 

145 

1.45006-03 

1.4450E-03 

l 

0.00100000 

0,84299956 

148 

1.48006-03 

1.4750E-03 

2 

0.00200000 

0.84499955 

149 

1.49006-03 

1.48506-03 

4 

0.00400000 

0.84899954 

151 

1.51006-03 

1.5050E-03 

3 

0.00300000 

0.85199954 

152 

1.52006-03 

1.51506-03 

3 

0.00300000 

0.85499954 

153 

1.53006-03 

1.52506-03 

l 

0.00100000 

0.85599954 

157 

1.57006-03 

1.56506-03 

l 

0.00100000 

0.85699953 

159 

1.59006-03 

1.58506-03 

3 

0.00300000 

0.85999953 

160 

1.60006-03 

1.59506-03 

l 

0.00100000 

0.86099952 

162 

1.62006-03 

1.61506-03 

1 

0.00100000 

0.86199952 

165 

1.65006-03 

1.64506-03 

1 

0.00100000 

0.86299951 

166 

1.66006-03 

1.65506-03 

2 

0.00200000 

0.86499951 

167 

1.67006-03 

1.6650E-03 

3 

0.00300000 

0.86799951 

169 

1.69006-03 

1.6850E-03 

2 

0.00200000 

0.86999951 

170 

1.70006-03 

1.69506-03 

2 

0.00200000 

0.87199950 

171 

1.71006-03 

1.7050E-03 

l 

0.00100000 

0.87299950 

176 

1.76006-03 

1.75506-03 

1 

0.00100000 

0.87399949 

100 

1.80006-03 

1.79506-03 

l 

0.00100000 

0.87499949 

181 

1.81006-03 

1.80506-03 

1 

0.00100000 

0.87599948 

183 

1.83006-03 

1.82506-03 

2 

0.00200000 

0.87799948 

184 

1.84006-03 

1.83506-03 

1 

0.00100000 

0.87899947 

185 

1.85006-03 

1.84506-03 

1 

0.00100000 

0.87999947 

186 

1.86006-03 

1.85506-03 

1 

0.00100000 

0.88099946 

188 

1.88006-03 

1.87506-03 

2 

0.00200000 

0.88299946 

39 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

189 

1.8900E-03 

1.8850E-03 

190 

1.9000E-03 

1.8950E-03 

191 

1.9100E-03 

1.9050E-03 

192 

1.9200E-03 

1.9150E-03 

194 

1.9400E-03 

1.9350E-03 

195 

1 .9500E-03 

1.9450E-03 

196 

1.9600E-03 

1 » 9550E-03 

198 

1.9800E-03 

1.9750E-03 

199 

1.9900E-03 

1.9850E-03 

200 

2.0000E-03 

l . 9950E-03 

202 

2.0200E-03 

2. 01506-03 

205 

2.0500E-03 

2.0450E-03 

207 

2.0700E-03 

2.0650E-03 

208 

‘•2.Q800E-03 

2.0750E-03 

209 

2.0900E-03 

2.0850E-03 

211 

2.1100E-03 

2.1050E-03 

213 

2 .13006-03 

2 . 1250E-03 

214 

2.1400E-03 

2. 1350E-03 

215 

2 . 1500E-03 

2 . 1450E-03 

216 

2.1600E-03 

2. 1550E-03 

218 

2 • 1800E-03 

2 . 1750E-03 

226 

2.2600E-03 

2.2550E-03 

227 

2.2700E-03 

2.2650E-03 

235 

2.3500E-03 

2.34506-03 

240 

2.4000E-03 

2.39506-03 

241 

2.4100E-03 

2.4050E-03 

243 

2.4300E-03 

2.42506-03 

244 

2.4400E-03 

2.43506-03 

248 

2.48OOE-03 

2.4750E-03 

252 

2.5200E-03 

2.5150E-03 

254 

2.5400E-03 

2.53506-03 

256 

2.56006-03 

2.5550E-03 

257 

2.5700E-03 

2 . 5650E-03 

258 

2.5800E-03 

2.5750E-03 

265 

2.65006-03 

2.64506-03 

266 

2.6600E-03 

2.6550E-03 

272 

2.72006-03 

2.71506-03 

275 

2.75006-03 

2.74506-03 

283 

2.83006-03 

2.8250E-03 

284 

2.84006-03 

2.83506-03 

292 

2.92006-03 

2.91506-03 

296 

2.96006-03 

2.95506-03 

302 

3. 02006-03 

3.0150E-03 

309 

3.09006-03 

3.08506-03 

317 

3. 17006-03 

3.16506-03 

329 

3.2900E-03 

3.2850En03 

332 

3. 3200E-03 

3.3150E-03 

333 

3.33006-03 

3.32506-03 

340 

3.40006-03 

3. 3950E-03 

348 

3.48006-03 

3.4750E-03 

354 

3. 54006-03 

3.53506-03 

356 

3.56006-03 

3.55506-03 

357 

3.5700E-03 

3. 5650E-03 

361 

3.6100E-03 

3.6050E-03 

362 

3 .62006-03 

3.61506-03 

365 

3.65006-03 

3.64506-03 

FREQ 

DENSITY 

CUMULATIVE 

3 

0.00300000 

0.88599946 

1 

0.00100000 

0.88699945 

1 

0.00100000 

0.88799945 

1 

0.00100000 

0.88899944 

2 

0.00200000 

0.89099944 

2 

0.00200000 

0.89299943 

2 

0.00200000 

0.89499943 

1 

0.00100000 

0.89599942 

1 

0.00100000 

0.89699942 

3 

0.00300000 

0.89999942 

1 

0.00100000 

0.90099941 

1 

0.00100000 

0.90199941 

1 

0.00100000 

0.90299940 

1 

0.00100000 

0.90399940 

2 

0.00200000 

0.90599939 

1 

0.00100000 

0.90699939 

1 

0.00100000 

0.90799938 

1 

0.00100000 

0.90899938 

1 

0.00100000 

0.90999937 

1 

0.00100000 

0.91099937 

1 

0.00100000 

0.91199936 

1 

0.00100000 

0.91299935 

1 

0.00100000 

0.91399935 

l 

0.00100000 

0.91499934 

1 

0.00100000 

0.91599934 

1 

0.00100000 

0.91699933 

2 

0.00200000 

0.91899933 

1 

0.00100000 

0.91999932 

1 

0.00100000 

0.92099932 

3 

0.00300000 

0.92399932 

2 

0.00200000 

0.92599931 

1 

0.00100000 

0.92699931 

i 

0.00100000 

0.92799930 

l 

0.00100000 

0.92899930 

i 

0.00100000 

0.92999929 

l 

0.00100000 

0.93099929 

l 

0.00100000 

0.93199928 

l 

0.00100000 

0.93299928 

l 

0.00100000 

0.93399927 

l 

0.00100000 

0.93499926 

l 

0.00100000 

0.93599926 

l 

0.00100000 

0.93699925 

l 

0.00100000 

0.93799925 

l 

0.00100000 

0.93899924 

l 

0.00100000 

0.93999924 

l 

0.00100000 

0.94099923 

i 

0.00100000 

0.94199923 

l 

0.00100000 

0.94299922 

3 

0.00300000 

0.94599922 

1 

0.00100000 

0.94699921 

1 

0.00100000 

0.94799921 

1 

0.00100000 

0.94899920 

1 

0.00100000 

0.94999920 

1 

0.00100000 

0.95099919 

1 

0.00100000 

0.95199919 

1 

0.00100000 

0.95299918 

40 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

366 

3.6600E-03 

3. 6550E-03 

369 

3.6900E-03 

3.6850E-03 

370 

3.7000E-03 

3.6950E-03 

376 

3.76006-03 

3.75506-03 

379 

3.7900E-03 

3.78506-03 

385 

3.85006-03 

3. 8450E-03 

386 

3. 8600E-03 

3.85506-03 

393 

3.93006-03 

3.92506-03 

405 

4«0500E“03 

4.04506-03 

417 

4.1700E-03 

4. 1650E-03 

444 

4.4400E-03 

4. 4350E-03 

448 

4.48006-03 

4.47506-03 

453 

4. 5300E-03 

4.52506-03 

458 

4.5800E-03 

4.57506-03 

463 

4.6300E-03 

4.6250E-03 

478 

4.7800E-03 

4.7750E-03 

489 

4.8900E-03 

4.88506-03 

522 

5.2200E-03 

5.21506-03 

525 

5.2500E-03 

5.24506-03 

528 

5.28006-03 

5.27506-03 

529 

5.29006-03 

5.28506-03 

537 

5.37006-03 

5.36506-0 J 

575 

5.75006-03 

5.74506-03 

579 

5.79006-03 

5.7850E-03 

597 

5.97006-03 

5.96505-U3 

610 

6.10006-03 

6.09506  03 

633 

6. 3300E-03 

6.3250E-03 

644 

6.4400E-03 

6.43506-03 

668 

6.68006-03 

6.67506-03 

690 

6.90006-03 

6.89506-03 

702 

. 7.0200E-03 

7.01506-03 

716 

7.16006-03 

7.15506-03 

735 

7.35006-03 

7.34506-03 

750 

7.50006-03 

7.49506-03 

044 

8.44006-03 

8.43506-03 

870 

8.70006-03 

8.69506-03 

902 

9.02006-03 

9.01506*03 

932 

9.32006-03 

9.31506-03 

971 

9 • ?',  00E-03 

9.70506-03 

973 

9.73006-03 

9. 72506-03 

982 

o, 82006-03 

9.81506-03 

983 

9.83006-03 

9,82506-03 

1001 

FREQ 

D6NSITY 

CUMULATIVE 

1 

S 00100000 

0.95399918 

1 

0. 00100000 

0.95499917 

1 

0.00100000 

0.95599917 

1 

0.00100000 

0.95699916 

1 

0.00100000 

0.95799915 

1 

0.00100000 

0.95899915 

l 

0.00100000 

0.95999914 

1 

0.00100000 

0.96099914 

- 1 

0.00100000 

0.96199913 

1 

0. OOIOOOOO 

0.96299913 

1 

0.00100000 

0.96399912 

1 

0.00100000 

0.96499912 

1 

0.00100000 

0.96599911 

1 

0.00100000 

0.96699911 

1 

0. OOIOOOOO 

0.96799910 

1 

0.00100000 

0.96899910 

1 

0.00100000 

0.96999909 

1 

0.00100000 

0.97099908 

1 

0.00100000 

0.97199908 

1 

0.00100000 

0.97299907 

l 

0. OOIOOOOO 

0.97399907 

1 

0.00100000 

0.97499906 

1 

0.00100000 

0.97599906 

L 

0. OOIOOOOO 

0.97699905 

1 

0.00100000 

0.97799905 

1 

0.00100000 

0.97899904 

1 

0.00100000 

0.97999904 

1 

0. OOIOOOOO 

0.98099903 

1 

0. OOIOOOOO 

0.98199902 

1 

0.00100000 

0.98299902 

1 

0.00100000 

0.98399901 

l 

0. OOIOOOOO 

0.98499901 

l 

0.00100000 

0.98599900 

1 

0.00100000 

0.98699900 

1 

0.00100000 

0.98799899 

1 

0. OOIOOOOO 

0.98899899 

1 

0.00100000 

0.98999898 

1 

0.00100000 

0.99099898 

1 

0. OOIOOOOO 

0.99199897 

1 

0.00100000 

0.99299897 

1 

0.00100000 

0.99399896 

1 

0.00100000 

0.99499895 

5 

0.00500000 

0.99999895 

TL-STAR  STATISTICS 


HISTOGRAM  GENERATOR  OUTPUT 
SAMPLE  SAMP  MEAN  SAMP  STO  DEV 

1000  1.1214906  03  4.1340416  02 


HIST  MEAN 
1.1313606  03 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

FREQ 

0 

0. 

0 

20 

4.0000E 

02 

3.90006  02 

2 

21 

4.2000E 

02 

4.1000E  02 

I 

22 

4.4000E 

02 

4.30006  02 

2 

23 

4.6000E 

02 

4.50006  02 

3 

24 

4.8000E 

02 

4.70006  02 

6 

25 

5.0000E 

02 

4.9000E  02 

5 

26 

5.2000E 

02 

5.10006  02 

7 

27 

5.4000E 

02 

5.30006  02 

9 

28 

5.6000E 

02 

5.50006  02 

11 

29 

5.8000E 

02 

5.70006  02 

12 

30 

6.0000E 

02 

5.90006  02 

9 

31 

6.2000E 

02 

6.1000E  02 

8 

32 

6.4000E 

02 

6.3000E  02 

10 

33 

6.6000E 

02 

6.5000E  02 

9 

34 

6.8000E 

02 

6.7000E  02 

17 

35 

7.0000E 

02 

6.90006  02 

15 

36 

7.2000E 

02 

7.1000E  02 

15 

37 

7.4000E 

02 

7.30006  02 

17 

38 

7.6000E 

02 

7.50006  02 

21 

39 

7.8000E 

02 

7.70006  02 

21 

40 

8.00006 

02 

7.90006  02 

21 

41 

8.2000E 

02 

8.10006  02 

19 

42 

8.4000E  .02 

8.30006  02 

20 

43 

8.6000E 

02 

8.50006  02 

25 

44 

8.80006 

02 

8.70006  02 

23 

45 

9 .00006 

02 

8.9000E  02 

19 

46 

9.20006 

02 

9.10006  02 

22 

47 

9.4000E 

02 

9.3000E  02 

17 

48 

9.6000E 

02 

9.50006  02 

18 

49 

9.8000E 

02  . 

9.7000E  02 

23 

50 

U 00006 

03 

9.9000E  02 

21 

51 

l .02006 

03 

1.01006  03 

23 

52 

1.04006 

03 

1.0300E  03 

18 

53 

1.06006 

03 

1.05006  03 

25 

54 

1.08006 

03 

1.0700E  03 

20 

55 

l. 10006 

03 

1.0900E  03 

15 

56 

1.12006 

03 

1.11006  03 

22 

57 

1.14006 

03 

1.1300E  03 

12 

58 

l. 16006 

03 

1.15006  ‘03 

27 

59 

1.18006 

03 

1.17006  03 

15 

DENSITY 

0.00000000 

0.00200000 

0.00100000 

0.00200000 

0.00300000 

0.00600000 

0.00500000 

0.00700000 

0.00900000 

o.onooooo 

0.01200000 

0.00900000 

0.00800000 

0.01000000 

0.00900000 

0.01700000 

0.01500000 

0.01500000 

0.01700000 

0.02100000 

0.02100000 

0.02100000 

0.01900000 

0.02000000 

0.02500000 

0.02300000 

0.01900000 

0.02200000 

0.0170000Q 

0.01800000 

0.02300000 

0.02100000 

0.02300000 

0.01800000 

0.02500000 

0.02000000 

0.01500000 

0.02200000 

0.01200000 

0.02700000 

0.01500000 


HIST  STD  DEV 
4. 128561E  02 

CUMULATIVE 

0.00000000 

0.00200000 

0.00300000 

0.00500000 

0.00800000 

0.01400000 

0.01900000 

0.02600000 

0.03500000 

0.04600000 

0.05800000 

0.06700000 

0.07500000 

0.08500000 

0.09400000 

0.11100000 

0.12600000 

0.14100000 

0.15799999 

0.17899999 

0.19999999 

0.22099999 

0.23999999 

0.25999999 

0.28499999 

0.30799999 

0.32699998 

0.34899998 

0.36599998 

0.38399998 

0.40699998 

0.42799998 

0.45099998 

0.46899997 

0.49399997 

0.51399997 

0.52899996 

0.55099996 

0.56299996 

0.58999995 

0.60499994 


interval 

UPPER  LIMIT 

CLASS  MARK 

FREQ 

DENSITY 

CUMULATIVE 

60 

1.2000E  03 

1.19006  03 

25 

0.02500000 

0.62999994 

61 

1.22006  03 

1.2100E  03 

12 

0.01200000 

0.64199994 

62 

1.24006  03 

1.2300E  03 

24 

0.02400000 

0.66599993 

63 

1.2600E  03 

1.2500E  03 

13 

0.01300000 

0.67899993 

64 

1.28006  03 

1.2700E  03 

13 

0.01300000 

0.69199993 

65 

1.3000E  03 

1.2900E  03 

14 

0.01400000 

0.70599993 

66 

1.3200E  03 

1 .31006  03 

17 

0.01700000 

0.72299992 

67 

1.34006  03 

1.3300E  03 

11 

0.01100000 

0.73399992 

68 

1.3600E  03 

1.3500E  03 

16 

0.01600000 

0.74999992 

69 

1.3800E  03 

1.3700E  03 

- 5 

0.00500000 

0.75499991 

70 

1.400QE  03 

1.3900E  03 

7 

0.00700000 

0.76199991 

71 

1.4200E  03 

1.4100E  03 

12 

0.01200000 

0.77399991 

72 

1.44006  03 

1.4300E  03 

11 

0.01100000 

0.78499991 

73 

1.4600E  03 

l .45006  03 

10 

0.01000000 

0.79499990 

74 

1.4800E  03 

1.47006  03 

9 

0.00900000 

0.80399990 

75 

1.50006  03 

1.49006  03 

10 

0.01000000 

0.81399990 

76 

1.5200E  03 

1.51006  03 

15 

0.01500000 

0.82899989 

77 

1.5400E  03 

1.53006  03 

3 

0.00300000 

0.83199989 

78 

1.5600F  03 

1.5500E  03 

11 

0.0110000G 

0.84299989 

79 

1.58006  03 

1.57006  03 

13 

0.01300000 

0.85599989 

80 

1.60006  03 

1.59006  03 

9 

0.00900000 

0.86499988 

81 

1.6200F  03 

1.6100E  03 

9 

0.00900000 

0.87399988 

82 

1.6400E  03 

1.63006  03 

13 

0.01300000 

0.88699988 

83 

1.66006  03 

1.6500E  03 

14 

0.01400000 

0.90099987 

84  • 

. 1.6800E  03 

1.67006  03 

10 

0.01000000 

0.91099987 

85 

1.70006  03 

1.69006  03 

3 

0.00300000 

0.91399987 

86 

1.72006  03 

1.71006  03 

6 

0.00600000 

0.91999987 

87 

1.74006  03 

1.73006  03 

9 

0.00900000 

0.92899986 

88 

1.76Q0E  03 

1.75006  03 

3 

0.00300000 

0,93199986 

89 

1.7800E  03 

1.77006  03 

3 

0.00300000 

0.93499986 

90 

1.8000E  03 

1.79006  03 

3 

0.00300000 

0.93799986 

91 

1.8200E  03 

1.81006  03 

2 

0.00200000 

0.93999986 

92 

1.8400E  03 

1.83006  03 

3 

0.00300000 

0.94299985 

93 

1.8600E  03 

1.85006  03 

4 

0.00400000 

0.94699985 

94 

1.8800E  03 

1.87006  03 

9 

0.00900000 

0.95599984 

95 

1.90006  03 

1.89006  03 

4 

0.00400000 

0.95999984 

96 

l .92006  03 

1 .91006  03 

2 

0.00200000 

0.96199983 

97 

1.94006  03 

1.93006  03 

1 

0.00100000 

0.96299983 

99 

1.9800E  03 

1.97006  03 

4 

0.00400000 

0.96699982 

100 

2.0000E  0? 

1.99006  03 

2 

0.00200000 

0.96899982 

101 

2.0200E  03 

2.01006  03 

l 

0.00100000 

0.96999981 

103 

2.0600E  03 

2*05006  03 

5 

0.00500000 

0.97499981 

105 

2.10006  03 

2.09006  03 

l 

0.00100000 

0.97599980 

106 

2.12006  03 

2.11006  03 

2 

0.00200009 

0.97799980 

107 

2.1400E  03 

2.1300E  03 

1 

0.00100000 

0.97899979 

108 

2.16006  03 

2.15006  03 

1 

0.00100000 

0.97999979 

109 

2.18006  03 

2.1700E  03 

1 

0.00100000 

0.98099978 

no 

2.2000E  03 

2.19006  03 

l 

0.00100000 

0.98199978 

111 

2.22006  03 

2.21006  03 

2 

0.00200000 

0.98399977 

112 

2.24006  03 

2.23006  03 

1 

0.00100000 

0.98499977 

113 

2.26006  03 

2.25006  03 

1 

0.00100000 

0.98599976 

114 

2.28006  03 

2.27006  03 

l 

0.00100000 

0.98699976 

118 

2.3600E  03 

2.35006  03 

l 

0.00100000 

0.98799975 

119 

2.38006  03 

2.37006  03 

l 

0.00100000 

0.98899975 

120 

2 .40006  03 

2.39006  03 

l 

0.00100000 

0.98999974 

122 

2.44006  03 

2.43006  03 

l 

0.00100000 

0.99099974 

43 


INTERVAL  UPPER  LIMIT 

123  2.4600E  03 

124  2.4800E  03 

128  2.5600E  03 

141  2.8200E  03 

163  3.2600E  03 

177  3.5400E  03 

501 


CLASS  MARK  FREQ  DENSITY 

2.4500E  03  1 0.00100000 

2.4700E  03  3 0.00300000 

2.5500E  03  2 0.00200000 

2.8100E  03  1 0.00100000 

3.2500E  03  l 0.00100000 

3.5300E  03  1 0.00100000 

0 0.00000000 


CUMULATIVE 

0.99199973 

0.99499973 

0.99699973 

0.99799972 

0.99899971 

0.99999971 

0.99999971 


44 
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U-STAR  STATISTICS 

HISTOGRAM  GENERATOR  OUTPUT  - • - • 

SAMPLE  SAMP  MEAN  SAMP  STD  DEV  HIST  MEAN 

1000  1.9693756-03  1.7030656-03  1.634699E-03 


INTERVAL 

UPPER  LIMIT 

CLASL  MARK 

FREQ 

DENSITY 

0 

1.0000E-03  ' 

0 

0.00000000 

9 

1.18006-03 

1.17006-03 

128 

0.12800000 

10 

1.2000E-03 

1 .19006-03 

66 

0.06600000 

11 

1.22Q0F-03 

1.2100E-03 

50 

0.05000000 

12 

1.2400E-03 

1.23006-03 

’ 40 

0.04000000 

13 

1.2600E-03 

1.25006-03 

47 

0,04700000 

14 

1.2800E-03 

1.27006-03 

29 

0.02900000 

15 

1.30008-03 

1.29006-03 

26 

0.02600000 

16 

1.32006-03 

1.3100E-03 

31 

0.03100000 

17 

1.34006-03 

1.33006-03 

23 

0.02300000 

18 

1.3600E-03 

1.35006-03 

23 

0.02300000 

19 

1.3800E-03 

1.37006-03 

25 

0.02500000 

20 

1.4000E-03 

1.39006-03 

22 

0.02200000 

21 

1.4200E-03 

1.41006-03 

22 

0.02200000 

22 

1.44G0E-03 

1.4300E-03 

8 

0.00800000 

23 

1.46006-03 

1.4500E-03 

16 

0.01600000 

24 

1.4800E-03 

i. 47006-03 

14 

0.01400000 

25 

1.50006-03 

1.49006-03 

8 

0.00800000 

26 

1.5200E-03 

1.51006-03 

15 

0.01500000 

27 

1.54006-03 

1.53006-03 

10 

0.01000000 

28 

1 . 5600E-03 

1.5500E-03 

18 

0.01800000 

29 

1.5800F-03 

1.57006-03 

10 

0.01000000 

30 

1.6000E-03 

1.59006-03 

7 

0.00700000 

31 

1.6200E-03 

1.61006-03 

7 

0.00700000 

32 

1 .64006-03 

1.63006-03 

10 

0.01000000 

33 

1.66006-03 

1.65006-03 

7 

0.00700000 

34 

1.68006-03 

1.67006-03 

10 

0.01000000 

35 

1.70006-03 

1.69006-03 

4 

0.00400000 

36 

1.72006-03 

1.7100E-03 

12 

0,01200000 

37 

1.74006-03 

1.73006-03 

2 

0.00200000 

38 

1.74006-03 

1.75006-03 

3 

0.00300000 

39 

1.78006-03 

1.77006-03 

4 

0.00400000 

40 

1.80006-03 

1.79006-03 

11 

0.01100000 

41 

1.82006-03 

1.81006-03 

4 

0.00400000 

42 

1 .84006-03 

1.83006-03 

4 

0.00400000 

43 

l. 86006-03 

1.85006-03 

7 

0.00700000 

44 

1.88006-03 

1.8700E-03 

9 

0.00900000 

45 

1.90006-03 

1.89006-03 

2 

0.00200000 

46 

1 .92006-03 

1.91006-03 

8 

0.00800000 

47 

1.94006-03 

1.93006-03 

5 

0.00500000 

48 

1. 96006-03 

1.95006-03 

4 

0.00400000 

HIST  STD  DEV 
8.2212586-04 

- tUMULATIVE 
0.00000000 
0.12800000 
0. 19400000 
0.24400000 
0.28399999 
0.33099999 
0.35999999 
0.38599999 
0.41699999 
0.43999999 
0.46299998 
0.48799998 
0.50999998 
0.53199998 
0.53999998 
0.55599997 
0.56999997 
0.57799996 
0.59299996 
0.60299996 
0.62099995 
0.63099995 
0.63799995 
0.64499995 
0.65499995 
0.66199995 
0.67199995 
0.67599994 
0.68799993 
0.68999993 
0.69299993 
0.69699992 
0.70799992 
0.71199992 
0.71599991 
0.72299991 
0.73199990 
0.73399990 
0.74199989 
0. 74699989 
0.75099988 


45 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

FREQ 

DENSITY 

CUMULATIVE 

49 

1.9800E-03 

1 . 9700E-03 

2 

0.00200000 

0.75299988 

50 

2.0000E-03 

1 .99006-03 

•6 

0.00600000 

0.75899988 

51 

2.0200E-03 

2.01006-03 

3 

0.00300000 

0.76199988 

52 

2.0400F-03 

2.0300E-03 

2 

0.00200000 

0.76399987 

53 

2.0600E-03 

2 . 0500E-03 

4 

0.00400000 

0.76799987 

54 

2.0800E-03 

2.0700E-03 

8 

0.00800000 

0.77599986 

55 

2 . 1000E-03 

2. 0900E-03 

3 

0.00300000 

0.77899986 

56 

2.1200E-03 

2. 1 100E-03 

3 

0.00300000 

0.78199986 

57 

2.1400E-03 

2. 1300E-03 

2 

0.00200000 

0.78399985 

58 

2.1600E-03 

2. 1500E-03 

2 

0.002C0000 

0.78599985 

59 

2.1800E-03 

2 • 1 700E-03 

6 

0.00600000 

0.79199985 

60 

2.2000E-03 

2. 1900E-03 

3 

0.00300000 

0. 79499^985 

61 

2.2200E-03 

2.2100E-03 

2 

0.00200000 

0.79699984 

62 

2.2400E-03 

2.2300E-03 

3 

0.00300000 

0.79999984 

63 

2.2600E-03 

2.25006-03 

2 

0.00200000 

0.80199984 

64 

2.2800E-03 

2.2700E-03 

3 

0.00300000 

0.80499984 

65 

2.3000E-03 

2.2900E-03 

4 

0.00400000 

0.80899983 

66 

2 . 3200E-03 

2.31006-03 

3 

0.00300000 

0.81199983 

67 

2.3400E-03 

2 • 3300E-03 

4 

0.00400000 

0.81599982 

69 

2.3800E-03 

2.37006-03 

9 

0.00900000 

0.82499982 

70 

2.4000E-03 

2.39006-03 

1 

0.00100000 

0.82599981 

71 

2. 4200E-03 

2.4100E-03 

2 

0.00200000 

0.82799981 

73 

2.46Q0E-03 

2.45006-03 

1 

0.00100000 

0.82899980 

74 

2.4800E-03 

2.47006-03 

l 

0.00100000 

0.82999980 

75 

2.5000E-03 

2.4900E-03 

2 

0.00200000 

0.83199979 

76 

2.5200E-03 

2.51006-03 

3 

0.00300000 

0.83499979 

77 

2.5400E-03 

2.53006-03 

2 

0.00200000 

0.83699979 

78 

2.5600E-03 

2.55006-03 

2 

0.00200000 

0.83899979 

79 

2.5800E-03 

2.57006-03 

2 

0.00200000 

0.84099978 

80 

2. 60006-03 

2.59006-03 

l 

0.00100000 

0.34199978 

81 

2.6200E-03 

2. 6100E-03 

l 

0.00100000 

0.84299977 

82 

2.6400E-03 

2.63006-03 

2 

0.00200000 

0.84499977 

83 

2.6600E-03 

2.65006-03 

4 

0.00400000 

0.84899976 

84 

2. 68006-03 

2.67006-03 

5 

0.00500000 

0.85399976 

85 

2.7000E-03 

2.69006-03 

2 

0.00200000 

0.85599975 

87 

2.74006-03 

2.73006-03 

1 

0.00100000 

0.85699975 

88 

2.7600E-03 

2.75006-03 

4 

0.00400000 

0.8609 **974 

90 

2.8000F-03 

2.7900E-03 

l 

0.00100000 

0.86199474 

91 

2.82006-03 

2.8100E-03 

3 

0.00300000 

0.36499973 

92 

2.84006-03 

2.8300E-03 

3 

0.00300000 

0.86799973 

93 

2. 86006-03 

2. 8500E-03 

2 

0.00200000 

0.  869^9973 

94 

2.88006-03 

2.87006-03 

3 

0.00300000 

87299973 

97 

2.94006-03 

2.93006-03 

l 

0.00100000 

(5. 87399972 

99 

2.98006-03 

2.97006-03 

2 

0.00200000 

0.87599972 

100 

3.00006-03 

2.99006-03 

3 

0.0030000P 

0.87899972 

101 

3.02006-03 

3.01006-03 

2 

0.0020000V 

0.880999^1 

102 

3.04006-03 

3.03006-03 

2 

0.00200000 

0-88299971 

103 

3.06006-03 

3.05006-03 

4 

0.00400000 

0.88699970 

104 

3.08006-03 

3.07006-03 

2 

0.00200000 

0.88099970 

105 

3.10006-03 

3.09006-03 

2 

0.00200000 

0.89099970 

106 

3.12006-03 

3.11006-03 

4 

0.00400000 

0. 69499969 

107 

3.14006-03 

3.1 300E-03 

1 

0.00100000 

0.89599968 

108 

3.16006-03 

3.1 500F-03 

4 

0.00400000 

0.89999968 

109 

3.18006-03 

3. 1700E-03 

1 

0.00100000 

0.90099967 

111 

3.22006-03 

3.21006-03 

1 

0.00100000 

0.90199967 

112 

3.24006-03 

3.23006-03 

2 

0.00200000 

0.90399966 
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INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

113 

3.2600E-03 

3.2500E-03 

114 

3.2800E-03 

3.2700E-03 

115 

3.3000E-03 

3.29006-03 

116 

3.32006-03 

3.31006-03 

117 

3.3400F-03 

3.33006-03 

121 

3.42006-03 

3.41006-03 

122 

3.4400E-03 

3.43006-03 

126 

3.52006-03 

3. 5100E-03 

128 

3.5600E-03 

3.55006-03 

129 

3.5800E-03 

3. 5700E-03 

130 

3.6000E-03 

3.59006-03 

132 

3.6400E-03 

3.63006-03 

134 

3.68006-03 

3.6700E-03 

135 

3. 70006-03 

3. 69006-03 

136 

3.72006-03 

3.71006-03 

137 

3.74006-03 

3.73006-03 

141 

3.82006-03 

3.8100E-03 

144 

3.88006-03 

3.87006-03 

146 

3.92006-03 

3.91006-03 

160 

4.00006-03 

3.99006-03 

154 

4.08006-03 

4.07006-03 

156 

4.12006-03 

4.  U006-03 

159 

4.18006-03 

4.17006-03 

163 

4.26006-03 

4.25006-03 

167 

4.34006-03 

4.33006-03 

173 

4. 4600E-03 

4.45006-03 

174 

4.4800F-03 

4.4700E-03 

175 

4,50006-03 

4,49006-03 

17S 

4.56006-03 

4.55006-03 

182 

4.64006-03 

4.63006-03 

105 

4.70006-03 

4.69006-03 

186 

4.72006-03 

4.71006-03 

187 

4.74006-03 

4.73006-03 

109 

4.78006-03 

4.77006-03 

191 

4.82006-03 

4.81006-03 

193 

4.86006-03 

4.85006-03 

197 

4.94006-03 

4.93006-03 

198 

4.96006-03 

4.95006-03 

200 

5.00006-03 

4.99006-03 

201 

FREQ 

DENSITY 

CUMULATIVE 

2 

0.00200000 

0.90599966 

1 

0.00100000 

0.90699966 

1 

0.00100000 

0.90799965 

2 

0.00200000 

0.90999965 

2 

0.00200000 

0.91199964 

1 

0.00100000 

0.91299964 

1 

0.00100000 

0.91399963 

1 

0.00100000 

0.91499963 

1 

0.00100000 

0.91599962 

* 1 

0.00100000 

0.91699962 

3 

0.00300000 

0.91999961 

l 

0.00100000 

0.92099961 

2 

0.00200000 

0.92299961 

2 

0.00200000 

0.92499960 

2 

0.00200000 

0.92699960 

2 

0.00200000 

0.92899960 

2 

0.00200000 

0.93099959 

1 

0.00100000 

0.93199959 

1 

0.00100000 

0.93299958 

2 

0.00200000 

0.93499958 

1 

0.00100000 

0.93599957 

l 

0.00100000 

0.93699957 

1 

0.00100000 

0.93799956 

1 

0.00100000 

0.93899956 

1 

0.00100000 

0.93999955 

1 

Q.OOIOQOOQ 

0.  94099955 

l 

0.00100000 

0,94199954 

1 

0.00100000 

0.94299953 

3 

0.00300000 

0.94599953 

1 

0.00100000 

0.94699953 

l 

G*  00100000 

0.94799952 

l 

0.00100000 

0.94899952 

1 

0.00100000 

0.94999951 

2 

0.00200000 

0.95199951 

2 

0.00200000 

0.95399950 

2 

0.00200000 

0.95599950 

1 

0.00100000 

0.95699950 

1 

0.00100000 

0.95799949 

2 

0.00200000 

0.95999949 

40 

0.04000000 

0.99999949 
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Z-ST  AT  I ST  ICS 


histogram  generator  output 


SAMPLE 

SAMP 

MEAN  SAMP  STD 

OEV 

HIST  MEAN 

HIST  STD  DEV 

1000  -3. 

206379E  00  1. 

0248866-01 

-3.1874296  00 

2. 253406E-01 

interval 

0 

UPPER  LIMIT 
-3.5000E  00 

CLASS  MARK 

FREQ 

4 

DENSITY 

0.00400000 

CUMULATIVE 

0.00400000 

1 

-3.4900E 

00 

-3.4950E 

00 

1 

0.00100000 

0.00500000 

3 

-3.4700E 

00 

-3.4750E 

00 

1 

0.00100000 

0.00600000 

4 

-3.46006 

00 

-3.4650E 

00 

1 

0.00100000 

0.00700000 

5 

-3.45006 

00 

-3.4550E 

00 

3 

0.00300000 

0.01000000 

i 6 

-3.4400E 

00 

-3.44506 

00 

7 

0.00700000 

0.01700000 

7 

-3.4300E 

00 

-3.43506 

00 

2 

0.00200000 

0.01900000 

8 

-3.4200E 

00 

-3.4250E 

00 

3 

0.00300000 

0.02200000 

9 

-3.4100E 

00 

-3.4150E 

00 

6 

0.00600000 

0.02800000 

10 

-3.40006 

00 

-3.4050E 

00 

5 

0.00500000 

0.03300000 

11 

-3. 3900E 

00 

-3.3950E 

00 

6 

0.00600000 

0.03900000 

12 

-3.38006 

00 

-3.3850E 

00 

9 

0.00900000 

0.04800000 

13 

-3.3700E 

00 

-3.3750E 

00 

10 

0.01000000 

0.05800000 

14 

-3.36006 

QO 

-3.3650E 

00 

8 

0.00800000 

0.06600000 

15 

-3.35006 

00 

-3.3550E 

00 

14 

0.01400000 

0.08000000 

16 

-3.34006 

00 

-3.3450E 

00 

15 

0.01500000 

0.09500000 

17 

-3.3300E 

00 

-3.3350E 

00 

13 

0.01300000 

0.10800000 

13 

-3.32006 

00 

-3.3250E 

00 

13 

0.01300000 

0.12100000 

19 

-3.3100E 

00 

-3.3150E 

00 

24 

0.02400000 

0. 14499999 

20 

-3*  3000E 

00 

-3.3050E 

00 

18 

0.01800000 

0.16299999 

21 

-3.29006 

00 

-3.29506 

00 

35 

0.03500000 

0.19799999 

22 

-3.28006 

00 

-3.28506 

00 

22 

0.02200000 

C. 21999999 

23 

-3.27306 

00 

-3.27506 

00 

22 

0.02200000 

0.24199999 

24 

-3.26006 

00 

-3.26506 

00 

31 

0.03100000 

0.27299999 

25 

-3.25006 

00 

-3.25506 

00 

28 

0.02800000 

0.30099999 

26 

-3.24006 

00 

-3.24506 

00 

44 

0.04400000 

0.34499999 

27 

-3.23006 

00 

-3.2350E 

00 

33 

0.03300000 

0.37T99999 

20 

-3.22006 

00 

-3.22506 

00 

34 

0.03400000 

0.41199999 

29 

-3.21006 

00 

-3.21506 

00 

30 

0.03000000 

0.44199998 

30 

-3.2000E 

00 

-3.20506 

00 

40 

0. 04000000 

0.48199990 

31 

-3.1900E 

00 

-3.19506 

00 

56 

0.05600000 

0.53799998 

32 

-3.18006 

00 

-3.18506 

00 

31 

0.03100000 

0.56899998 

33 

-3.17006 

00 

-3.17506 

00 

43 

0.04300000 

0.61199997 

34 

-3.16006 

00 

-3.16506 

00 

40 

0.04000000 

0.65199997 

35 

-3.15006 

00 

-3.15506 

00 

36 

0.03600000 

0.68799997 

36 

-3.14006 

00 

-3.14506 

00 

46 

0.04600000 

0.73399997 

37 

-3.1300E 

00 

-3.13506 

00 

24 

0.02400000 

0.75799996 

30 

-3.12006 

00 

-3.12506 

00 

22 

0.02200000 

0.77999996 

39 

-3.11006 

00 

-3.11506 

00 

18 

0.01800000 

0. 79799996 

40 

. -3.10006 

00 

-3.10506 

00 

35 

0.03500000 

0.83299996 

41 

-3.0900E 

00 

-3.09506 

00 

19 

0.01900000 

0.85199995 

/ 


INTERVAL 

UPPER  LIMIT 

CLASS  MARK 

FREQ 

DENSITY 

CUMULATIVE 

42 

-3.0800E 

00 

-3.0850E 

00 

31 

0.03100000 

0.88299995 

43 

-3.0700E 

00 

-3.0750E 

00 

19 

0.01900000 

0.90199994 

44 

-3.0600E 

00 

-3.0650E 

00 

20 

0.02000000 

0.92199994 

45 

-3.0500E 

00 

-3.05506 

00 

14 

0.01400000 

0.93599994 

46 

-3.0400E 

00 

-3.0450E 

00 

13 

0.01300000 

0.94899993 

47 

-3.0300E 

00 

-3.0350E 

00 

11 

0.01100000 

0.95 99 999 3 

48 

-3.0200E 

00 

-3.02506 

00 

6 

0.00600000 

0.96599993 

49 

-3.0100F 

09 

-3.0150E 

00 

12 

0.01200000 

0.97799993 

50 

-3.000QE 

00 

-3.00506 

00 

9 

0.00900000 

0.98699992 

51 

-2.9900E 

00 

-2.9950E 

00 

' 4 

0.00400000 

0.99099991 

52 

-2.9800E 

00 

-2.9850E 

00 

4 

0.00400000 

0.99499991 

54 

-2.9600E 

00 

-2.9650E 

00 

2 

0.00200000 

0.99699990 

55 

-2.9500E 

00 

-2.9550E 

00 

2 

0.00200000 

0. 99899990 

57 

-2.9300E 

00 

-2.9350E 

00 

1 

0.00100000 

0.99999990 

101 

0 

0.00000000 

0.99999990 
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APPENDIX 


The  primary  data  on  tube  life  which  forms  the  basis  of  the 
assertion  that  F(t)  is  log-normal  was  obtained  for  the  M113  cannon. 

This  is  given  in  table  2.-.  The  24  failure  times  (cycles  of  max  rated 
pressure)  are  listed  by  rank  together  with  the  median  rank  statistic 
for  a sample  of  24. 

The  median  ranks  were  plotted  versus  failure  times  on  log-normal 
graph  paper  and  on  Weibull  graph  paper  in  graphs  8 and  9,  respectively. 

As  suggested  by  the  appearance  of  the  data  on  log-normal  graph  paper, 
a conservative  estimate  of  the  underlying  distribution  was  drown  thru 
the  data.  Points  from  the  lower  tail  of  estimated  distribution  were 
plotted  on  Weibull  paper  with  the  experimental  data  to  provide  additional 
visual  indication  of  tha  tenablllty  of  the  log-normality  assumption 
vercua  sltsrnatlve  assumptions.  For  example,  the  estimated  points  seem 
to  fair  more  smoothly  with  tha  axparimantal  points  than  would  a straight 
llna  on  Weibull  paper. 
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TABLE  1 


Failure  Times 

Median  Rank 

1 

2194 

.0289 

2 , 

2196 

.0698 

3 

2826 

.1108 

4 

3025 

.1518 

5 

3169 

.1927 

6 

3309 

.2337 

7 

3323 

.2747 

8 

3438 

.3156 

9 

3449 

.3566 

10 

3476 

.3976 

11 

3606 

.4385 

12 

3698 

.4795 

13 

3847 

.5205 

14 

4451 

.5615 

15 

4671 

.6024 

16 

4685 

.6434 

17 

4697 

.6844 

18 

4766 

.7253 

19 

5035 

.7663 

20 

5397  . 

.8073 

21 

6192 

.8482 

22 

7238 

.8892 

23 

7631 

.9302 

24 

9652 

.9711 

eve* 

4272.68 
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